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We study a non-Hermitian generalization of quantum systems in which an imaginary 
vector potential is added to the momentum operator. In the tight-binding approximation, 
we make the hopping energy asymmetric in the Hermitian Hamiltonian. In a previous article, 
we conjectured that the non-Hermitian critical point where the energy gap vanishes is equal 
to the inverse correlation length of the Hermitian system and we confirmed the conjecture 
for two exactly solvable systems. In this article, we present more evidence for the conjecture. 
We also argue the basis of our conjecture by noting the dispersion relation of the elementary 
excitation. 
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1. Introduction 

In this paper, we study a non-Hermitian generalization of strongly correlated quantum 
systems in which an imaginary vector potential ig (where ^ is a real vector) is added to the 
momentum operator. The non-Hermitian kinetic energy in the continuous space is given by^ 

W. = t^±^. (1, 

Im 

Its second-quantized form within the tight-binding approximation is given by 

= E E {^'-'^'^i+e.^^ + . (2) 

In this article, we focus on the one-dimensional case with a constant imaginary vector potential 
and hereafter use 

= -tE(^'4+iCx + e-54c,+i), (3) 



where 5 is a real constant; in other words, we make the hopping energy asymmetric. More 
je] 

CxCx+n by e~"'^ in the original Hermitian Hamiltonian. 



generally, we multiply the right hopping energy cl._^_^Cx by e"^ and the left hopping energy 
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The purpose of the non-Hermitian generahzation is to obtain a length scale inherent in 
the wave function of the Hermitian system only from the non-Hermitian energy spectrum. 
The non-Hermitian generalization was first applied to the one-electron Anderson model by 
Hatano and Nelson.^ Their model is, in one dimension, 

L L 
Wrandom(5) =t ^ (e» |x + 1) + c'^ + 1|) + ^ |x) (x| , (4) 

x=l x=l 

where is a random potential at site x and we require the periodic boundary condition. As 
we increase the non-Hermiticity g, a pair of eigenvalues collide at a point g = gc and then 
become complex. It was revealed^ that the non-Hermitian critical point gc is equal to the 
inverse localization length of the eigenfunction of the original Hermitian Hamiltonian. 

We here apply the non-Hermitian generalization to systems without randomness but 
with interactions. We conjectured in the previous article^ that we can obtain the correla- 
tion length from the non-Hermitian generalization of strongly correlated quantum systems; 
the non-Hermitian critical point g = gc where the energy gap from the ground state collapses 
is equal to the inverse correlation length of the ground state of the Hermitian system. We con- 
firmed the conjecture for two one-dimensional exactly solvable systems: the Hubbard model 
in the half- filled case and the S = 1/2 antiferromagnetic XX Z chain in the Ising-like region. 

In the present article, we give further evidence for the conjecture in various levels of 
certainty. We also clarify the reason why we can obtain the inverse correlation length by our 
non-Hermitian generalization, by noting the dispersion relation of the elementary excitation. 
We show that the non-Hermitian generalization is actually equivalent to replacing k with k+ig 
in the dispersion relation of the elementary excitation; we thus seek a zero of the dispersion 
relation by computing the non-Hermitian critical point gc- 

In § 2, we confirm our conjecture for the following exactly solvable systems in one dimen- 
sion: the S = 1/2 ferromagnetic isotropic XY chain in a magnetic field in § 2.1; the half-filled 
Hubbard model in § 2.2; the S" = 1/2 antiferromagnetic XXZ chain in the Ising-like region in 
§ 2.3; the Majumdar-Ghosh model in § 2.4. 

In § 3, we numerically analyze non-Hermitian models of finite size L. We calculate the 
non-Hermitian "critical" point gc{L) where the energy of the eigenstate corresponding to the 
ground state in the limit L ^ oo becomes complex; we then obtain an extrapolated estimate 
g'c(oo). We numerically confirm that the estimate gdoo) and the inverse correlation length of 
the Hermitian systems are consistent for the Hubbard model and for the XXZ model. We also 
analyze an unsolved model, namely S = 1/2 antiferromagnetic Heisenberg chain with nearest- 
and next-nearest-neighbor interactions. 

In the summary, we conclude our discussions by giving a remark on applying our non- 
Hermitian generalization. 
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2. Non-Hermitian analysis of solvable models 

For four exactly solvable one-dimensional systems, we here confirm our conjecture that the 
non-Hermitian critical point is equal to the inverse correlation length of the Hermitian system. 
We also reveal in this section that the non-Hermitian generalization results in replacing the 
real momentum k with a complex one k + ig in the dispersion relation of the elementary 
excitation. That is, the non-Hermitian generalization makes the Hermitian Hamiltonian of 
the form 

^= E ^^^'^ivk (5) 

—1T<k<1T 

transformed to 

n{9)= E <k + ig)rilrik. (6) 

— 7r<fc<7r 

Finding the non-Hermitian critical point where the energy gap above the ground state collapses 
is thereby equivalent to seeking a zero of the dispersion relation in the complex momentum 
plane. The zero of the dispersion relation may be related to the inverse correlation length 
of the Hermitian system; see Appendix B. We confirm the equivalence for the isotropic XY 
chain in § 2.1 and for systems solved by the Bethe-ansatz method in § 2.2 and § 2.3. We 
thereby show for the three models that the non-Hcrmitian critical point gc where the energy 
gap vanishes is rigorously equal to the inverse correlation length. 

We also analyze in § 2.4 the Majumdar-Ghosh model, for which we do not know the energy 
gap exactly; only approximate estimates are known. We still show that the non-Hermitian 
critical point gc where the approximate estimates of the energy gap vanish are equal to the 
inverse correlation length calculated by finite-size scaling of the correlation function of the 
ground state of the Hermitian Majumdar-Ghosh model. Wc also argue that the non-Hermitian 
critical point gc obtained from the approximate estimates is the exact one. 

2.1 The isotropic XY chain in a magnetic field 

We first consider the non-Hermitian generalization of the S = 1/2 ferromagnetic isotropic 
XY chain in a magnetic field. The Hermitian Hamiltonian of this model is 

Hxy = ~jY1 i^iSf+i + sfsi,) - h Si, (7) 

1=1 1=1 
where we set J > 0. The Hamiltonian (7) is transformed into 

= X^(4+iCi + cfci+i) + ^ - ^ 4^1 (8) 

1=1 1=1 
by the Jordan- Wigner transformation 

cj = {-2y-'StSI . . . S]_,S-, c] = (-2)^-iSf 5| . . . (9) 
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(a) (b) 

Fig. 1. The eigenvalue distributions of the non-Hermitian isotropic XY chain iov h> J for (a) g = 
and (b) < g < Qc- The symbol x indicates the ground state. 



Applying the non-Hermitian generalization of the form (3) to the Hamiltonian (8), we have 

nxYig) = -^J2(^'cl^ici + e-^cjci+i) + — -hJ2cW (10) 
;=i 1=1 

By the inverse Jordan- Wigner transformation, the Hamiltonian (10) is transformed back into 

L L 



1=1 1=1 

= -jj^ [cosh5 (5f5f+i + 5f+J - isinh^ (Sf^f+i + SfSl,)] - hj^ Sf. (11) 

1=1 1=1 

We can immediately diagonalize the non-Hermitian Hamiltonian (10) in the momentum 
space in the form 

'HxY{g)= J2 <k + ig)clck + ^ (12) 

— 7r<fc<7r 

with the Fourier transformation 

c. = ^E^-''-" ct = -Lfje^«4, (13) 

where €{k) = — J cos k — h. The non-Hermitian generalization thus shifts the momentum k by 
ig in the dispersion relation of the elementary excitation. 

We define the non-Hcrmitian critical point gc as the point where the energy gap above the 
ground state vanishes. For h > J, all cigcnvahics are real at the Hcrmitian point g = and 
e(0) gives a finite energy gap (Fig 1. (a)). As wc turn on the non-Hermiticity g, all eigenvalues 
except for k = and ibvr immediately spread into the complex k plane (Fig. 1 (b)). The energy 
gap AE(g) in the limit L — > oo is given by the elementary excitation energy in the form 

AE{g) = e{±TT + ig) = h- Jcoshg. (14) 
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Fig. 2. The non-Hermitian critical point Qc of the isotropic XY chain. The ground-state critical point 
of the Hermitian system is h/J= 1. 




The energy gap vanishes at 




which is the non-Hermitian critical point of the model. We stress here that eq. (15) gives a 
zero of the dispersion relation: e(ib7r + igc) = 0. It is plausible that the zero of the dispersion 
relation in the complex k plane gives a typical length scale of the system. Indeed, we have 

e(±7r + V6=0. (16) 

Figure 2 shows the non-Hermitian critical point Qc as a function of h/J. Hermitian system 
is gapless (the XY phase) for h < J and hence we have the non-Hermitian critical point Qc = 
in the region. We can immediately confirm that the analytical expression of the non-Hermitian 
critical point gc is equal to the inverse correlation length of the Hermitian system, which was 
obtained by the quantum transfer matrix method.^ 

2.2 The half-filled Hubbard chain 

We next consider the non-Hermitian generalization of the one-dimensional Hubbard model 
in the half-filled case: 

L L 
Wcharge(5) = -^Y. {'''^ 4+1 ,a^i + e-^C^^q+i,^) + ^ Yl ^^^^ 

1=1 o-=T,J, 1=1 
where U is the on-site Coulomb repulsive interaction. Fukui and Kawakami^ solved the non- 
Hermitian model (17) by the Bethe-ansatz method. They showed that we eliminate the charge 
gap, namely the Hubbard gap, as we increase the non-Hermiticity g. The fact is understood 
intuitively by computing the first-order perturbation of the Hamiltonian (17) with respect to 
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9- 

^charge (ff) ^charge (0) 

L 

= -9t Y,(4^i,fi,^ - + cj^i^^q,^ - = -ig{J^ + J;), (18) 

1=1 

where Ja{(^ =T) i) is the paramagnetic current operator defined by 

L 

= -itY,i4+l,a^l,'7 - 4,a^l+l,<7)- (19) 

1=1 

Hence J| + is the charge current operator; see eq. (28) below for comparison. At the 
Hermitian point g = oi the Hamiltonian (17) in the limit L oo, the ground state is the 
Mott insulator for any finite U. As we increase the non-Hermiticity g, that is, as wc try to 
produce the charge fiow (19) in one direction, wc expect that the ground state changes from 
the Mott insulator to an extended state, or a "metallic" state at a point g = gc\ this is the 
non-Hermitian critical point where the Hubbard gap vanishes. 

In order to solve the non-Hcrniitian Hamiltonian (17), we make the ansatz for the right 
eigenfunction 5'^, which is the same as in the Hermitian case: 

N 

*^(a;i,a;2, • • .,xm\xm+i, ■ ■ ■,xn) = ^[Q, P] exp(z ^ fcp^XQ^), (20) 

P j=i 

where M and N are the number of the down spins and the number of the electrons, respec- 
tively. The Bethe-ansatz wave function (20) is different from that of Fukui and Kawakami's; 
their Bethe-ansatz wave function is given by replacing kj with kj—ig in eq. (20). The difference 
lies only in the definition of kj . We here use the form (20) in order to compare the dispersion 
relation (21) below with the Hermitian one. We put the down spins at xi,X2, ■ ■ ■ ,xm and 
the up spins at xm+i, ■ ■ , xn- The symbols P = [Pi, Pi-, ... , Pn) and Q = {Qi, Q2, ... , Qn) 
are two permutations of the set (1,2, . . . ,N) with 1 < xq^ < xq^ < • • • < xq,^ < L. The 
symbol [Q, P] is a set of N\ x N\ coefficients depending on the two permutations P and Q. 
The "dispersion relation" with respect to the quasimomentum kj is given by 

e{kj) = -2tcos{kj + ig), (21) 

where kj is a solution of the Bethe-ansatz equations: 

^ sm{kj+ig)-Af^ + iU/At 
'''^^'''^^^=} Linik,+^g)-A,-^U/4t (^ = 1'---'^)' 

^ smjkj + ig) - Ag + iU/4t _ ^ - Af^ - iU/At 

f}^si^{kj + ig)-A^-iU/At ll^A^-Ap + iU/At ^ ^ ^ ) 

The non-Hermitian generalization thus results in replacing kj with kj + ig in the dispersion 
relation. 
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Fig. 3. The g dependence of the Hubbard gap A for U/t = 2. The non-Hermitian critical point is 
5c = 0.038. 



Fukui and Kawakami'^ obtained an analytical expression of the non-Hermitian critical 
point gc in the limit L ^ oo, in the form 

f°° X + iU/At 

gc = arcsinh(?7/4t) + 2i arctan — — — — a{X)dX, (23) 

J —oo U / 'it 

where the distribution function cr(A) of the spin rapidity A is given by 

I rco fU \ 
a{X) = — / sech — uj cos{Xoj)Jo{uj)(kj (24) 
27r Jq \4i J 

with Jq{uj) the Bessel function of the first kind. After some algebra in Appendix A, we can 
show that the non-Hermitian critical point gc is actuahy equal to the inverse correlation length 
1/^ due to the charge excitation of the Hermitian Hubbard model in the form 

^ - arcsinh - 2 H M^) ^^^H{U ^ .25) 

which was obtained by Stafford and Millis^ from finite-size scaling of the Drude weight. 

We note here that eq. (25) is equal to the imaginary part of the momentum at a zero of 
the dispersion relation of the charge excitation in the complex momentum space. This is the 
same situation as in (16) for the XY chain; see Appendix B.l for details. 

Figure 3 exemplifies how the Hubbard gap collapses as we increase the non-Hermiticity g. 
The ground-state energy does not change and the Hubbard gap gradually decreases before it 
vanishes at g = gc- The way the energy gap collapses is different from that for the Anderson 
model discussed by Hatano and Nelson,^ where the energy gap decreases almost suddenly as 
g gets close to g = gc- Unfortunately, it is difficult to know the ground-state properties of 
the non-Hermitian Hubbard model (17) for g > gc- However, we expect that the ground-state 
energy becomes complex in the region g > gc on the basis of finite-size data shown in § 3.1. 
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2.3 The 5 = 1/2 antiferromagnetic XXZ chain 

We next develop a parallel discussion for a spin system, namely the S = 1/2 antiferro- 
magnetic XXZ chain. The non-Hermitian Hamiltonian of the XXZ chain is given by^ 

L 



=1 L 



(26) 



for J > 0. We set S^^^ = hereafter. The non-Hermitian Hamiltonian Hxxzig) of 
A = 1 is derived as an effective Hamiltonian of the following non-Hermitian Hubbard model: 

L L 
1=1 1=1 

(27) 

We have the Hamiltonian (26) with A = 1 by considering the second-order perturbation with 
respect to the kinetic energy around the degenerate ground state. Note the difference between 
eqs. (17) and (27); the first-order perturbation with respect to the non-Hermiticity g gives 

Wspm(<7) - :^spin(0) = -ig{J^ - Ji), (28) 

where J| — is the spin current operator. Thus we expect that the non-Hermiticity g induces 
a spin current and eliminates the spin gap. 

Albertini et al.^ solved eq. (26) in the massive region A > 1 in the subspace S^^^. = 0. 
In order to solve the non-Hermitian Hamiltonian (26), we make the following ansatz for the 
right eigenfunction which is the same as in the Hermitian case: 

M 

^'^(xi,X2,...,xm) = ^[Q,-P]exp(z^A;p^XQ^), (29) 

P j=i 

where we put the down spins at xi,X2, ■ ■ ■ ,xm {1 < xi < X2 < ■ ■ ■ < xm < L). The Bethe- 
ansatz wave function (29) is different from that of Albertini et aZ.'s; their Bethe-ansatz wave 
function is given by replacing k with A; — 2ig in eq. (29); see a comment below eq. (20). 
The "dispersion relation" with respect to the quasimomentum kj is given by 

e{kj) = J [A + cos(% + 2ig)] , (30) 

where kj is a solution of the Bethe-ansatz equations: 

exM L) - (-1)^-' TT exp[z(fc,- + k) - 4ff] + 1 - 2A exp(zfe,- - 2g) . . _ ^ 2 • • • M) 
exp(zfc,L)-( 1) [[^,,^^i^kj + k)-4g] + l-2AeMikl-29) ' 

(31) 

The non-Hermitian generalization results in replacing kj by kj + 2ig in the dispersion relation. 

Albertini et al. obtained an analytical expression in the limit L — > 00 of the non-Hermitian 
critical point gc at which the spin gap vanishes, in the form^ 

^ _ 7 , (-l)"tanh(n7) 
5c - 2 + 2^ n ' 

n=l 
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where 7 = arccoshA. Although Albcrtini et al. did not point out the fact, the expression (32) is 
actuaUy weU known as the inverse correlation length'^ at the Hermitian point 5 = obtained by 
the quantum transfer matrix method. We also note that eq. (32) is related to the imaginary 
part of the momentum at a zero of the dispersion relation of the spinon excitation in the 
complex momentum space; sec Appendix B.2 for details. 

It is again difficult to know the ground-state properties of the non-Hermitian XX Z chain 
(26) for g > Qc- However, wc expect that the ground-state energy becomes complex in the 
region g > on the basis of finite-size data shown in § 3.2. 



2.4 The Majumdar- Ghosh model 

We now consider the non-Hermitian generalization of the Majumdar-Ghosh model^'^ 

Hug = jY^^liStiSf + S+Sf_^,) + SfSf^,] + ^ ^^[^(S+^Sr + S+Sf_^,) + SfSf^,]. 
1=1 1=1 

The Hamiltonian (33) has two- fold degenerate ground states and has a finite energy gap.^° 
However, only variational estimates of the energy gap are known. -"^^'^^ We introduce the non- 
Hermitian Hamiltonian of the Majumdar-Ghosh model by referring to the non-Hermitian 
generalization of the XXZ chain discussed in § 2.3: 

^ 1 

1=1 



■ =1 



(33) 



The above Hamiltonian can be derived from the effective Hamiltonian of the non-Hermitian 
t-t'-U model as follows: 

L 

1=1 
L 

1=1 

L 

+ t/^4^Q,t4;Q,i. (34) 
1=1 

Let us first calculate the correlation length of the Hermitian Majumdar-Ghosh model 
from finite-size scaling of the correlation function. The two-fold degenerate ground states of 
the Hermitian Majumdar-Ghosh model (33) are 

= ; ^ =(^I + ^11), 

V2 + 4 • 2-V2 
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where wave functions $i and $n are defined by 

$1 = 01,2 O 4>3,A (8) • • • (8) 4)L-l,L, 

^11 = 02,3 <8) 04,5 <8) • • • <8) 0L,1 (36) 

with 0j^j denoting the singlet state of a pair of spins at sites i and j: 

= ^(1 Uj) - I iiW)- 

Note that $i and $i are not orthogonal: ($i|$ii) = ($ii|$i) = 2"^/^+^ The states (35), on 
the other hand, are orthonormal. 

The correlation functions {SqS^)± = {'^±\SqS^\'^±) with respect to the two-fold degen- 
erate ground states and of the system of size L are given by 

— 1 m O-L/2+2 1 , 

iSoSlh = -ii^rj^ = ~T 2-^/^ + O ((2-/2)2 
(SoSkh = ^ ^ = ±2-V2 + o ((2-^/2)2) (for i > 1), 

(^o^2Vi)± = ^ J2-L/2+1 = ^2-^/' + O ((2-^2)2^ (fo^ . > (37) 
Assuming finite-size correction of the correlation function in the form 

{S^S',)l = {S^S^)oo + O (exp(-L/0) , (38) 
we obtain the correlation length 

We next calculate the non-Hermitian critical point gc where approximate estimates of the 
energy gap vanishes. 

The dispersion relation was first obtained by Shastry and Sutherland^2 with a trial wave 
function and then by Caspcrs et a/.^^ with a variational wave function. The dispersion relation 
of the Hermitian Majumdar-Ghosh model given by Shastry and Sutherland^2 ig 

€s{k) = J (^^ + cosk^ ; (40) 

the excitation energy at A; = ivr determines the energy gap AE^- When we turn on the 
non-Hermiticity g as in eq. (33), the momentum k in eq. (40) is replaced by A: + 2ig. The 
dependence of the energy gap AEs{g) on the non-Hermiticity g is 

AEsig) = es(±7r + 2ig) = j(^l- cosh(25)^ . (41) 

As shown in Fig. 4, the non-Hermitian critical point g^ where the energy gap AEg vanishes 
is gc = In 2/2. 
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The variational form of the dispersion relation given by Caspers et al. is 
(60r + 34) cos 2k + (139r + 170) cos k - (475r - 731) 
^ ~ (8r + 4) cos 2k - (28r + 70) cos k - (52r - 86) ^ ' 



with 

v'2cos 2/c-20cosA; + 43 



(43) 



5 + 4 cos k 

the excitation energy at A; = ±Tr determines the energy gap AEq- When we turn on the 

non-Hermiticity g, the momentum k in eq. (42) is again replaced by A; + 2ig. The dependence 

of the energy gap AEc{g) on the non-Hermiticity g is 

{60R + 31) co^h(-lg) - (139i? + 170) cosh(2g) - ( i75i? - 731) 
^^^^ (8i? + 4) cosh(45) + (28i2 + 40) cosh(25) - (52i? - 86) ^ ' 

with 



/ 2cosh(4g) + 20cosh(2g)+43 
Y (5-4cosh(25))2 • ^^""^ 

As shown in Fig. 4, the non-Hermitian critical point gc is gc = In 2/2 again. In both cases, the 
non-Hermitian critical point gc is equal to the inverse correlation length = In 2/2. 

We now give a piece of evidence that the non-Hermitian critical point is exactly = In 2/2 
by showing that the ground-state property dramatically changes at g = In 2/2. The non- 
Hermitian Majumdar- Ghosh model (33) has two-fold degenerate dimer ground states as in 
the Hermitian case. The right eigenfunctions and the left eigenfunctions (^'^| of the 

ground states of the system of size L are given by 
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where the wave functions l^^i^n) and jj| are defined by 





= 10^) ^ 


3 10^4) <^ • ■ 


■■^\<1>l-i,l), 




= 10^3) <^ 


5 10^5) ® • ■ 






= {4'l2\ ^ 








= ('/'2,3l ^ 




••<^('/'L,ll 



(47) 

with Ict^fj) and ((?!)^'^| denoting weighted singlet states of a pair of spins at sites i and j: 



|<^5.) = i=(e-^|Td,)-e^|i.T,)) 
1 



{^lj = —{e^Ui,\-e-^U,\). (48) 

Equation (47) is consistent with eq. (36) transformed by a many-body version of the imaginary 
gauge transformation^ 

iIj{xi,X2,...;9) =e<^^i'''ilj{xi,X2,...;0). (49) 

Note that = (^>fi|${}) = 1 and = ($h|$f ) = {e^^ + e"^^) /2^/^ but 

and (^'±1 are bi-orthonormal. 

The correlation functions of the non-Hermitian system, {SqS^)± = {^\\SqS^\"^^) with 
respect to the two-fold degenerate ground states \"^^) and are obtained in the forms 

e9L + g-flL 



-1t2- 



2^/2 



^^"^ 2V2 



2i/± 



-gL 



2^/2 



=9^ 



2±2- 



+ e- 



-gL 



(for i > 1), 



2^/2 



+ 



(^g^2Vi)±= (fori>l). (50) 

2 + 2 — 

^ ^ 2^/2 

Figure 5 shows the correlation function in the limit L — ^ 00 in the region 0<5<ln2/2 and in 
the region 5 > In 2/2, respectively. The ground state is dimerized in the region 0<5<ln2/2 
and is an extended state in the region g > In 2/2. The phase transition from the dimer 
state to the extended state reminds us of the localization-delocalization transition of the non- 
Hermitian random Anderson model discussed by Hatano and Nelson.^ The phase transition 
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point gc = In 2/2 is then naturally regarded as the non-Hermitian critical point. We thus 
confirm from the above discussions that the non-Hermitian critical point is equal to the 
inverse correlation length of the Hermitian systems for the Majumdar- Ghosh model. 

3. Numerical analysis of non-Hermitian models 

In the previous section, we discussed the non-Hermitian generalization of exactly solvable 
models and we confirmed the conjecture that the non-Hermitian critical point Qc where the 
energy gap vanishes is equal to the inverse correlation length of the Hermitian system. We 
now numerically show that the inverse correlation length is consistent with the extrapolated 
estimate gdoo) of finite-size data gc{L) where an eigenvalue which corresponds to the ground 
state in the limit L ^ oo, becomes complex. We show the above for the Hubbard model in § 3.1, 
for the S = 1/2 antiferromagnetic XXZ model in § 3.2 and for a frustrated quantum spin chain 
in § 3.3. Although we do not know the correlation length of the S = 1/2 antiferromagnetic 
Heisenberg chain with nearest- and next-nearest-neighbor interactions, we show in § 3.3 that 
the numerical estimate gdoo) is consistent with the ground-state phase diagram. 
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Fig. 6. (a) The spectrum flow of the eigenvalues per site around the ground state for L = 4 with 
U/t = 2 as we increase the non-Hermiticity g which ehminates the charge gap. (b) The 1/L plot 
ofgdL). 



3.1 The non-Hermitian Hubbard chain 

We first analyze the non-Hermitian Hubbard models (17) and (27) of size L. We define 
the non-Hermitian "critical" point gc{L) of a finite system as the point where the energy 
gap between the ground state and a low-lying excited state vanishes and beyond which the 
ground-state energy becomes complex. We here show that the extrapolated estimate gdoo) 
of finite-size data gc{L) is close to the exact value of the correlation length. 

We first use the non-Hermitian generalization of the form (17) in the subspace S^^^ = 0, 
which eliminates the charge gap. All eigenvalues are real at the Hermitian point 5 = 0. Upon 
increasing g, a pair of eigenvalues move on the real axis. They spread into the complex plane 
when g exceeds a value gc{L)- 

Figure 6 (a) shows the spectrum flow of the eigenvalues per site for L = 4 around the 
ground state for U = 2t. The eigenvalues of the ground state and the third excited state 
move toward each other on the real axis and spread into the complex plane as soon as the 
two eigenvalues collide, which gives the non-Hermitian "critical" point gc{L)- The eigenvalues 
of the first and the second excited states scarcely move. The movement of the ground-state 
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energy is presumably a finite-size effect; the ground-state energy does not change for g < Qc 
for the infinite system as shown in Fig. 3. 

We numerically estimated gc{L) for L = 4,6,8 and 10 and extrapolated them to gc{oo) 
by considering the finite-size correction as follows: 

g^{L) = 5e(oo) + a/L + 0(1/L2). (51) 

Figure 6 (b) shows the 1/L plot of {L) ; this implies that we have to consider different finite- 
size corrections in the case L = 4n (for L = 4 and 8) and in the case L = 4n -|- 2 (for L = 6 
and 10). In Fig. 6 (b), the line A is the linear fitting of gc{L) for L = 6 and 10 and the line B 
is that for L = 4 and 8; both lines are determined by the least-squares method under the 
condition that they have the same intercept gc{oo). The final estimate of gc{oo) is 0.037, while 
the Bethe-ansatz method yields = VCchaxge — 0.038. Our estimate is consistent with the 
exact value. It is quite remarkable to obtain such a good estimate from data for such small L. 

Figure 7 (a) shows the spectrum fiow for L = 4 around the ground state when we use the 
non-Hermitian generalization of the form (27), which eliminates the spin gap. The eigenvalues 
of the first and the second excited states move toward each other, while the eigenvalues of 
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the ground state and the third excited state scarcely move. Wc presume that the energy gap 
between the ground state and the third excited state is caused by the charge excitation while 
the energy gap between the first and the second excited states is caused by the spinon exci- 
tation. We expect that the ground state and the first excited state are eventually degenerate 
in the thermodynamic limit. Hence we regard the collision of the first and second excited 
states as the ground-state transition. This behavior implies the charge-spin separation of one- 
dimensional quantum systems in the low-energy region. Figure 7 (b) shows the 1/L plot of 
gc{L). We obtain the extrapolated estimate gdoo) by the same least-squares method as we 
used above; the line A shows the linear fitting of gc{L) for L = 4 and 8 and the line B is that 
for L = 6 and 10, which yields the extrapolated estimate gdoo) = —0.003. Our estimate is 
also consistent with the exact value gc = 1/Cspinon = 0. 

3.2 The non-Hermitian XX Z chain 

Next, we numerically calculate the non-Hermitian "critical" point gc{L) of the XX Z chain 
of size L. We obtain the extrapolated estimate gdoo) of finite-size data gd^) and show that 
the estimate gd^o) is consistent with the inverse correlation length of the Hermitian system 
as was for the Hubbard model. Figure 8 (a) shows the spectrum flow per site around the 
ground state of the XXZ chain of L = 8 as we increase the non-Hermiticity ^ for A = 3. 
The pair of the first and the second excited states undergoes the real-complex transition, 
which is the same as in Fig. 7 (a). Note that the ground state and the first excited state of 
Hermitian finite systems are eventually degenerate in the thermodynamic limit. We therefore 
expect that the real-complex transition point of the first and second excited states converges 
to the non-Hermitian critical point gc in the limit L ^ oo. We thereby use the real-complex 
transition point in order to define the non-Hermitian "critical" point gdL). 

We then extrapolate the finite-size data. Figure 8 (b) shows the extrapolation of gd^) for 
the XXZ chain with A = 3. The extrapolated estimate gd^o) by linear fitting in the form 



for L = 12, 14 and 16, is 0.231. In order to take the finite-size data gd^) foi^ small L into 
consideration, we also calculate the extrapolated estimate gd^o) by a second-order fitting in 
the form 



for L = 4, 6, . . . , 14, 16 to obtain gdoo) = 0.235. Both estimates are consistent with the inverse 
correlation length g^ = 1/^ = 0.244 calculated analytically. 

3.3 The NNN Heisenberg chain 

We consider the non-Hermitian generalization of the S = 1/2 antiferromagnetic Heisenberg 
chain with nearest- and next-nearest-neighbor interactions; we hereafter call this model NNN 



gdL) =gd^) + a/L + 0(1/1^) 



(52) 



gdL) = gdoo) + a/L + t/L^ + 0(1/l3) 



(53) 
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Fig. 8. (a) The spectrum flow of the eigenvalues per site of the XXZ model for L = 8 with A = 3. 
(b) The finite-size plot of gc{L).^ 



Heisenberg chain. The Hermitian Hamiltonian of this model is 

L 

Hnnn = J ^[^1 ■ + aSi ■ S/+2], (54) 
1=1 

where J > and a > 0. We require the periodic boundary conditions. At the point a = 0, 
the model is the standard Heisenberg chain. The ground state is a spin fluid state and the 
energy gap is zero. At the point a = 1/2, the model is the Majumdar- Ghosh model^'^ and 
the energy gap is flnite. Okamoto and Nomura^^ numerically showed that a massive-massless 
transition occurs at Oc = 0.2411. 

We here analyze the non-Hermitian Hamiltonian of the NNN Heisenberg chain 

L L 

1=1 1=1 

(55) 

We numerically estimated the non-Hermitian "critical" point gc{L) of finite systems in the 
subspace S^^^ = where an eigenvalue which corresponds to the ground state in the ther- 
modynamic limit becomes complex. Figure 9 shows the spectrum flow of the eigenvalues per 
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Fig. 9. The spectrum flow of the real part of the eigenvalues per site around the ground state for 
L = 8 with (a) a = 0.49, (b) a = 0.5 and (c) a = 0.51. 



site around the ground state for L = 8 with a = 0.49, a = 0.5 (the Majumdar-Ghosh point) 
and a = 0.51. In Fig. 9 (a) for a = 0.49, as we increase the non-Hermiticity g, the energy 
gap between the first excited state, which corresponds to one of the degenerate ground states 
in L 00, and the second excited state, which corresponds to the first excited state in 
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L oo, vanishes at g = gci{L) = 0.24. These two eigenvalues become complex in the region 
gci{L) < g < gc2{L) before they become real again at g = gc2{L) = 0.42. We define the 
non-Hermitian "critical" point of a finite system for a < 0.5 as the point g = gci{L) where 
the first excited-state energy first becomes complex. 

In Fig. 9 (b) for a = 0.5, the two-fold degenerate ground states exist for any g and the 
energy gap between the ground state and the first excited state vanishes at g = gc{L) = 0.35. 
However, these eigenvalues do not become complex for g > gc{L). This is presumably because 
gci{L) = gc2{L) = gc{L). We still regard gci{L) as the non-Hermitian "critical" point gdL) 
for a = 0.5. 

In Fig. 9 (c) for a = 0.51, the ground state and the first excited state never become 
complex for any g. The ground state for a > 0.5 is an incommensurate state of the spiral 
phase. Our non-Hermitian generalization of the form (55) is presumably not appropriate 
for detecting the incommensurate state in the region a > 0.5; we may need other types of 
non-Hermitian generalization. 

We hereafter restrict ourselves to the region < a < 0.5. We extrapolate the finite-size 
data gdL) for L = 4, 8, 12 and 16 by the linear fitting in the form 

gdL) =gd^) + a/L + 0(1/ L^) (56) 

and by the second-order fitting in the form 

gdL) = 5c(oo) + a/L + b/L^ + 0(1/^3). (57) 

Figure 10 shows the extrapolated estimates yc(oo) in the region < a < 0.5. The second-order 
estimate </c(oo) around a = is almost zero and is consistent with = in the limit L oo. 
The extrapolated estimates gdoc) at a = 0.5 are around 0.35 for both linear and second-order 
fitting, which is consistent with the inverse correlation length = In 2/2 (= 0.347); see 
§. 2.4. The estimate gdoo) is almost zero in the region < a < 0.25 and is finite in the region 
0.25 ^ a < 0.5. This is consistent with the massive-massless transition^*^ at a = 0.2411 if we 
admit that 5c(oo) is equal to the inverse correlation length. It is remarkable that gdco) is at 
least approximately equal to the inverse correlation length all through the region < a < 0.5 
though we do not know the exact dispersion relation of the elementary excitation of this 
model. 

4. Summary and discussions 

To summarize, we considered the non-Hermitian generalization of strongly correlated 
quantum systems, where we make the hopping energy of the original Hermitian Hamiltonian 
asymmetric. We rigorously showed that the non-Hermitian critical point where the energy gap 
vanishes is equal to the inverse correlation length of the Hermitian systems for the ferromag- 
netic isotropic XY chain in the magnetic field, the half-filled Hubbard model, the 5 = 1/2 
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Fig. 10. The extrapolated estimates gdoo) as a function of a by (a) a linear fitting and by (b) a 
second-order fitting. 



antiferromagnetic XXZ chain in the Ising-hke region and the Majumdar- Ghosh model. We 
argued that the non-Hermitian generalization replaces the real momentum k with the com- 
plex one k + ig or k + 2ig in the dispersion relation of the elementary excitation. By the 
non-Hermitian generalization, we effectively seek the zero of the dispersion relation in the 
complex momentum space; the imaginary part of the momentum at the zero is related to the 
inverse correlation length of the Hermitian system. 

We also discussed the 5=1/2 antiferromagnetic Heisenberg chain with the nearest- and 
the next-nearest-neighbor interactions, which is unsolved analytically. We presented numeri- 
cal evidence that the extrapolated estimates of the non-Hermitian "critical" point for finite 
systems is consistent with the inverse correlation length. 

The non-Hermitian generalization of the form (3) is not always appropriate for having 
the inverse correlation length. The model (55) with a > 0.5 is an example. The 5 = 1/2 
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ferromagnetic transverse Ising chain 

H = -jY,SfSf^,-hY,Sf (58) 

1=1 1=1 

is another example. The elementary excitation of eq. (58) is given by 

?7fc = cos OkCk - sin ^fecl^., -q-k = sin 6'fec|. + cos 0kC-k, (59) 

where 



9k = — arctan 



J sin A; 



J cos k + 2h 



(60) 



It has the dispersion relation 



€{k) = ^{Jcosk/2 + + (Jsinfc/2)2. (61) 

By replacing k with k + ig in the dispersion relation e(/c), we obtain the non-Hermitian 
Hamiltonian in the form (6). This is transformed back to the spin Hamiltonian of the form 

oo L 

= Yl E(-2)'^^f+l • • • SUn-l [iPn - rn)Sr^nSr + {Pn + rn)SlJ', + ^gn(5f+,5f - Sf^^Sf 
n=l 1=1 

-P^j^St, (62) 

where coefficients Pn,Qn and are given by the following integrals: 

1 r fJcosk A e{k + ig) + e{k - ig) , 

= ^ £ eik + ig)-eik-^g) ^^^^^^^^^ 



J sink e{k + ig) + e{k — ig) 



s\n{nk)dk. (63) 



This non-Hermitian Hamiltonian is very complicated; interactions among spins beyond the 
nearest neighbor sites emerge as soon as g is finite. It is because its elementary excitation is 
obtained by the Bogoliubov transformation; the creation and annihilation operators at two 
different wave numbers k and —k are mixed. Conversely, the non-Hermitian generalization of 
the simple form (3) does not produce the dispersion relation e{k + ig) in this model. We may 
need another principle of non-Hermitian generalization. 
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Appendix A: Equality of and 1/^ for the Hubbard model 

Wc show for the Hubbard model that the non-Hermitian critical point in eq. (23) and 
the inverse correlation length 1/^ in eq. (25) are actually equal. The non-Hermitian critical 
point Qc is given by 



9c = lim 

A^oo 



arcsin 



inh(C//4i) + '^^ j arctan ^^^^^^j^a{\)d\ 



= lim 

A— »oo 



arcsinr 



5mh{U/4:t) + — d\ arctan 
71" J -A 

Using the variable transformation 



U/4:t 

X + iU/At cos(cjA)Jo(u;) 



U/4:t Jo cosmU/4:t)uj] 
e = arctan(A/(i7/4t) + i) 



duo 



(A-l) 



(A-2) 



with 



we have 



5c = lim 

A— ►oo 



= lim 

A— voo 



tan0i = -^+z, tan02 = ^ + i, 

cos((C//4i)w tan 61 - z(J7/4t)a;) Jo(u;)(iw 



. , , I /'^^ {UlAt)e 
arcsinh(?7/4t) + - ' ^ ' > 



arcsm. 



h(C//4t) + - / 



n {u/At)e r 

vr Jg^ cos2 9 [Jq 

{U/At)Jo{Lj)dLO 

cosh{{U/At)io) 



cosh((C//4t)w) 



(A-3) 



de 



sm{{U/4:t)uj tan 9 - i{U/At)oj) 
{UlAt)uj ' 



6*2 /■6»2 



Oi J til 



sin((?7/4t)a;tan6' - i{U/4:t)Lo) 



de 



(A-4) 



We rewrite the term Ii in the form 
1 



{U/U)u 
1 



[^2sin((J7/4t)u;tane2 - i{U/At)uj) - 0i sin((C//4t)a; tan^i - i{U/4t)u!)] 



h + 6*2) sin (u; A), 



{U/4t)uj 

where the coefficients 9i and $2 for A 3> 1 take the form 



with |(5i|, 1^21 ^ 1- Because of 

tan $1 

we have 



^l = -|-^l, 02 = ^-62 



11.11 

~ — , tan 02 



tan Si Si 



tan S2 S2 ' 



(A5) 



(A6) 



(A-7) 



1 



{U/At)uj 
1 



Si — S2) sin(a;A) 

U/At U/At 



{U/At)u V i{U/At) - A i{U/At) + A 



sin(a;A) 
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2i{U/4t) 



uj{{U/Uf + A2) 



sin(a;A) 



0. (A-8) 

Next we calculate the integral I2 in eq. (A-4). By using the variable transformation x = 
iaxiO — i, we have 

rA/{U/4t) sm{{U/4t)ux) dx 

I2 = lim / 

-A/(l//4t) 



A^'^ J-A/(u/4t) {U/U)u 1 + (a; + if 



f 

J — ( 



sin((;7/4t)a;x) _ vr 



dx = ^ (1 - e-2-(t^/4t)). (A.9) 



(?7/4t)u;(x2 + 2ix) 2i{U/it)u 
We thus arrive at 

5c =arcsinh(C//40 T -^^W^.^T^i^ ' e-'-^^l''^)d.^ 

=arcsmh(C//40 - 2 ^^^^^^J^^A^^d... (A-IO) 

We thus confirmed that the expression (A- 10), or eq. (23) is actually equal to the inverse 
correlation length of the charge excitation (25). 

Appendix B: A zero of the dispersion relation of the elementary excitation and 
the correlation length 

In § 2.1, we showed cq. (16) for the XY chain. In this section, we argue similar relations 
for the half-filled Hubbard model and the 6* = 1/2 antiferromagnetic XXZ chain. That is, 
the imaginary part of the momentum at a zero of the dispersion relation of the elementary 
excitation in the complex momentum space is related to the correlation length. 



B. 1 The half-filled Hubbard chain 

We consider the Hermitian Hubbard model 

L L 

i=i,(7=T,i 1=1 
in the half-filled case. The charge excitation at the quasimomentum k = kh has the excitation 
energy in the form 

£ih) = U + 4tcoskn + 8tJ^ -S-—^L^du^ (B.2) 

and has the momentum in the form 

p{k,) = k, + 2j^ ^^^^^u./2t^ d^. (B-3) 

By analytic continuation of the dispersion relation of the Hermitian model, we search a zero of 
£{kh) along the axis Kep{kh) = ±7r in the complex momentum space, which gives the Hubbard 
gap. Fukui and Kawakami^ pointed out that S{kh) vanishes at kh = =b7r-|-zarcsinh([/"/4t). The 
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p{kh) = ±Tr + i 
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At 



Jo{uj)smh{U/4t) 



doj 



(B-4) 



/o uo{l + e^u/2t^ 

We note that the imaginary part of the momentum (B-4) is equal to the inverse correlation 
length 1/^ given in eq. (25). That is, we have p{kh) = ±7r + i/^ when £{kh) = 0. This is the 
same situation as in eq. (16) for the XY chain. 

B.2 The S = 1/2 antiferromagnetic XX Z chain 

We next consider the S = 1/2 antiferromagnetic XXZ chain 



^ n 



(B-5) 



in the region A > 1. The spinon excitation at the quasimomenta A = Xg^ and A = Agj has the 
excitation energy in the form 

27r Jsinh7 



£{\s^-i A, 



S2J 



7 



(p(A.i) + p(A,,)), 



where the distribution function p(A) of the rapidity A is given by 



2 cosh(n7) 

and we set 7 = arccoshA. The momentum of the spinon excitation takes the form 

7(A^, +Xsi) . sin(n7Asi) + sin(n7A52) 
P{'^si,'^s2) — ^ <- 2^ — 



n=l 



n cosh(n7) 



(B-6) 



(B-7) 



(B-8) 



By analytic continuation of the dispersion relation of the Hermitian model, we search a zero 
of £{Xsi , Xs^) along the axis Rep(A5^ , A52) = iTr in the complex momentum space, which gives 
the spinon energy gap. Albertini et al.^ pointed out that £{Xsj^,Xs2) vanishes at A^^ = Xg^ = 
±7r/7 + i. The momentum p(Asi, A52) at this point is 

(—1)" tanh(n7) 



^(Asi.Asj) = ±Tr + i 



7 



n=l 



n 



(B-9) 



We note that the imaginary part of the momentum (B-9) is equal to twice the inverse correla- 
tion length 1/^ given in eq. (32). That is, we have j3(Asi, A^j) = ±7r-|-2i/^ when f (A^j , A^j) = 0. 
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